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Abstract 

Extending a previous single-temperature model, an electrostatic gyrofluid model that includes 
anisotropic temperatures (Tn ^ T±) and can treat general nonlinear situations is constructed. The 
model is based on a Lagrangian formulation of gyrofluid dynamics, which leads to an exact energy 
conservation law. Diamagnetic cancelations are inserted manually in such a way that energy 
conservation is preserved. Comparison with previous models shows a very good agreement for 
zero-Larmor-radius terms in the gyrofluid equations of motion. 



PACS numbers: 52.35.Ra, 52.30.Ex, 52.65.Tt 



I. INTRODUCTION 



The gyrofluid model [l]-[6] is one of the most commonly used tools for understanding and 
explaining tokamak turbulent transport. Although it is less accurate than the gyrokinetic 
approach, it is computationally more efficient and economic, and it retains most of the 
qualitative features necessary to explain turbulent transport [2] -[6]. 

To correctly compute turbulent transport, and especially for those computations that 
treat large-amplitude disturbances, or those that treat strong spatial variations in the plasma 
parameters, a gyrofluid model should satisfy an exact energy conservation law. In linear 
theory, a small deviation from energy conservation leads to a small change in the growth rate. 
Turbulence, however, also includes several conservative transfer pathways that influence the 
overall statistical equilibrium. Violation of energy conservation in this case can lead to 
a source of free energy, and, hence, large errors in the description of the turbulent state. 
Moreover, to study the nonlinear evolution subsequent to a linear growth rate, the model 
should be fully nonlinear, while conserving energy. 

In a previous paper [7], we derived a fully nonlinear set of electrostatic gyrofluid equa- 
tions that satisfy an exact energy theorem. This set contains the evolution equations for 
density, momentum, and perpendicular temperature, as well as the polarization equation 
for the electrostatic potential. The energy conservation law, which is guaranteed by the 
Noether method through the Lagrangian formulation of gyrofluid dynamics, is a unique and 
important feature of this model, which is thought to be especially important for the numer- 
ical simulation of turbulence. Although our previous single-temperature gyrofluid model is 
consistent, it does not contain the parallel-temperature effects necessary for a quantitatively 
accurate description of turbulence in magnetized plasmas. In the present paper, we extend 
our previous model to include these effects. Although the model is still electrostatic, it 
should give an accurate description of low (5 turbulent plasmas. 

The extension to anisotropic temperatures (Tj| ^ Tj_) turns out to be nontrivial. In the 
single-temperature case, the constraints used in the variation of the Lagrangian density can 
be chosen to be the particle and entropy conservation laws. In the two-temperature case, one 
needs an additional constraint and one can no longer rely on (only) conservation laws. The 
approach used here is to derive the constraints for the parallel and perpendicular pressure 
through the full Chew-Goldberger-Low (CGL) equations [8, 9]. A further difference with our 
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previous work is the use of the virtual fluid displacement £ in deriving Eulerian variations 
for the gyrofluid moments, instead of the Lagrangian-multiplier method (which turns out to 
be algebraically more involved). 

Gyrofluid equations that include anisotropic parallel and perpendicular pressures have 
been derived previously through the derivation of moments of the gyrokinetic equation [1]- 
[6]. Although these equations do not conserve energy, they agree to a large extent with our 
present gyrofluid model. 

The remainder of this paper is organized as follows. In Sec. II, we present the gyrofluid 
Lagrangian and give the constraints on the density and anisotropic pressures used in the 
variational principle. In Sec. Ill, we present the variation of the gyrofluid action functional, 
and derive the gyrofluid equations of motion before the insertion of diamagnetic cancela- 
tions, as well as the polarization equation. In Sec. IV, we derive the energy conservation law 
by applying the Noether method. In Sec. V, we describe how the diamagnetic cancelations 
are inserted into our gyrofluid equations and present the final equations for density, momen- 
tum, parallel and perpendicular pressures. We also present the explicit form of the energy 
conservation law with each energy-exchange term clearly identified. Lastly, in Sec. VI, we 
summarize our work and present our conclusions. 

II. GYROFLUID LAGRANGIAN DENSITY AND LAGRANGIAN 
CONSTRAINTS 

A. Gyrofluid Lagrangian Density 

The four-moment gyrofluid Lagrangian for the anisotropic-temperature model is con- 
structed in the same way as that of the one-temperature model [7], based on the work of 
Pfirsch and Correa-Restrepo [10]. Here, the gyrofluid Lagrangian density is defined as 

C f = \nm (u\ + |u £ | 2 ) - (p ± + |) + en (A." - ||0||) , (1) 

where u E = (c/B) b X V0 is the E x B velocity, u\\ = u ■ b denotes the gyrofluid velocity 
parallel to the magnetic field B = VxA = 5b, and the nonlocal operator ||...|| denotes a 
gyrofluid gyro-averaging operator defined as follows. First, we note that, in our four-moment 
model, the electrostatic gyro-potential ||0|| depends on the perpendicular temperature Tj_ = 
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p_i_/n and is assumed not to depend on higher-order gyrofluid moments. Thus, a general 
expression for ||0|| can be given as 

U\\ - v vl) = £ ^ p 2fe viV, (2) 

where p 2 = (mc 2 / e 2 B 2 )T±_ is the squared Larmor radius and V(x) represents a smooth 
function such that V(0) = 1, i.e., we have ||0|| — > in the zero-Larmor-radius (ZLR) limit. 
Two choices for V(x) will be discussed below; here, we note that the operator in the 
Taylor expansion (2) acts on alone. The Hermitian conjugate operator ||...||^ is defined in 
terms of Eq. (2) as 

where x is an arbitrary function and we made use of the identity 

xHW = Elixir + v ± - [(xi + vlx2 + - -^v^ - 0v ± (xi + vix 2 + •••) 

+ (X2 + •••)V±V 2 L - V10V ± (X2 + •••) + •••], (4) 

where 

v = g(xJMD ^ W (Q) v , fc 

Xfe ~ 9(Vf0) 2 fe A;! ' 
Note that the operator in Eq. (3) now acts on x an d powers of the Larmor radius p. 

The Ex B kinetic term m|u^| 2 /2 in Eq. (1) is the long-wavelength form of gyro-screening 
(i.e., the ZLR part of the second-order potential term in the gyrokinetic Hamiltonian [11]), 
which plays a major role in the derivation of the polarization equation. Although the E x B 
kinetic term is only an approximation of the full second-order gyrokinetic Hamiltonian, it 
is useful since the gyro-screening correction to the potential is important only for large 
flow velocities and long wavelengths. Thus, the addition of full finite-Larmor-radius (FLR) 
effects for this quadratic term will not alter any computational results. The fact that the 
gyro- screening term appears only in the Hamiltonian leads to the appearance of polarization 
effects only in Poisson's equation for the scalar potential 0, which is desirable for computa- 
tional reasons. 

Although the drift-fluid Lagrangian of Ref. [10] and our gyrofluid Lagrangian are similar, 
there exist three differences: (a) the gyrofluid Lagrangian (1) contains FLR corrections to 
the scalar potential 0, (b) the anisotropic temperatures T± and T| are treated separately 



rather than as a single isotropic temperature T, and (c) we have arranged the terms such 
that the gyro-screening term appears only in the Hamiltonian (i.e., not in the symplectic 
part involving contraction with u; see Ref. [7] for further details). The total Lagrangian 
density for our four-moment gyrofluid model is the sum of the gyrofluid Lagrangian density 
(1) and the Lagrangian density of the electric field, expressed as 

C = l -nm (u\ + \u E \ 2 ) - n(r ± + ^ + era (a • " - ||0||) + H (5) 

where the variational fields are the four gyrofluid moments (ra, u,p\\,p±) for each fluid species 
(here, summation over particle species is assumed wherever appropriate) and the electro- 
static potential 0. 



B. Lagrangian Constraints 

In order to proceed with our variational principle, we introduce constraints on the gy- 
rofluid moments (n,p\\,p±), based on the continuity (mass conservation) equation 

y +nV'U = fl, (6) 

at 

where d/dt = d/dt+u • V denotes the total time derivative, and the modified CGL equations 
for the perpendicular and parallel pressures [8, 9]: 

-P- + pii V-u + 2pn bb : Vu = - 2 V-qin + 4qin • (b • Vb), (7) 

dt 

^ + p ± V • u + p ± (I - bb) : Vu = - V • qj_± - 2 q,| ± • (b • Vb), (8) 

where qyj. and qj_j_ are the parallel and perpendicular heat fluxes, respectively. The pressure 
constraints, without the terms including the heat fluxes, can be also be obtained by general 
conservation laws. If the heat fluxes are removed from the pressure equations (7)- (8), all 
gyrofluid quantities rf = (n,p\\,p±) are directly related to the velocity u. Using these 
equations as constraints for the Lagrangian density (5), we can obtain the evolution equation 
for the gyrofluid velocity u. 

We obtain expressions for the Lagrangian variations Ar] a in terms of the virtual fluid 
displacement £ by taking the limits 

lim | ^- At] = An a and lim (u At) = £, 



so that the variations of the pressures from Eqs. (7) and (8) become 

Ap|| = V-£ - 2 P]l bb : V£ 



> . 



(9) 



Ap ± = -p ± V-£ - p± (I - bb) : V£ 

The heat-flux terms in Eqs. (7)-(8), which were not used in the Lagrangian variations (11) 
for p\\ and p±, are added to the pressure equations later on since they play an important role 
in the diamagnetic cancelations; these cancelations refer to the fact that the diamagnetic 
velocity does not lead to advection. The present Lagrangian formulation also requires one 
additional Lagrangian variation: 



Au = di/dt = d£/dt + u-V£ 



(10) 



for the fluid velocity u. For a correct treatment of FLR effects associated with the gyrofluid 
electrostatic potential ||0||, instead of the pressure variations we will use the temperature 
variations, such that the set of the Lagrangian variations used will be: 

An = -nV-£ 



Au = dt/dt + u-VZ 



ATI i 



AT, 



- 2Tji bb : V£ 



=- -Ti (I - bb) : V£ 



(11) 



where An is obtained from Eq. (6). 

Lastly, one should note that the gyrofluid Lagrangian density (5) was constructed, not 
derived, and that the method adopted here is not a mathematical derivation of the gyrofluid 
Lagrangian from the gyrokinetic (single-particle) Lagrangian. The validity of the present 
gyrofluid Lagrangian follows from the validity of the resulting gyrofluid equations, as well 
as the energy conservation law. The same applies to the constraint equations, e.g., Eqs. (7)- 
(8), which can be introduced arbitrarily, and their validity follows from the validity of the 
resulting evolution equations. 
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III. GYROFLUID DYNAMICAL EQUATIONS 



In this Section, we present a variational principle based on the action functional S = 
/ Cd 3 xdt using the virtual-displacement method associated with the Eulerian variations 
(5n,5u,5T\\,5T±). Here, each Eulerian variation 5x is defined in terms of its Lagrangian 
variation Ax as 8\ = A\ — £ ■ V%, so that, using Eq. (11), we find 

5n = - V • (n £) 



5u = d t £ + (u-V)£ - (£-V)u 



> . 



(12) 



«JT|| = -| • VT|, - 2T|, bb : V£ 

5T ± = - V • (T ± £) + Tj_ bb : V£ 

From the Eulerian variation of the gyrofluid Lagrangian (5), we can now derive the moment 
equation for the gyrofluid velocity u, as well as the polarization equation for the scalar 
potential 0. Note that the Eulerian variation 5T± in Eq. (12) is used also in connection with 
the FLR dependence of 



6\\<f>\\(<f>,T ± ) ee + ST ± 



9HW 

dT ± 



where 



dT ± 



2e 2 B 2 \ 2 x 



V10 = fi(0,T ± ;5), 



(13) 



(14) 



and V'(x) denotes the first derivative of V(x). As a result of the mass scaling in Eq. (14), 
the ion FLR correction f2j is much larger than the electron FLR correction Q e . 



A. Variation of the Lagrangian Density 

It is straightforward to derive the total variation 8C of the gyrofluid Lagrangian density 
(5), which depends on the variational fields (n, u, Tji, T±, 0), so that we obtain 



/ dC \ dC dC r n 

5C = (W0 - enll^llj + 5n — + Su-— - [-57], + (1 + eQ) nST ± 



(15) 
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Here, from Eq. (5), we find 



dC 

dn 



m 
~2 



n ( m My b + - A ) , (16) 



where we introduced the definition 



e^T ±] B) = eU\\ - (mc 2 /2 B 2 ) | V0| 2 , 



(17) 



which combines the linear (with FLR corrections) and nonlinear (in the ZLR limit) electro- 
static terms of the gyrokinetic Hamiltonian [11]. 

By inserting the Eulerian variations (12) into the gyroffuid Lagrangian variation (15), we 
obtain (after some algebra) 

5C _ _ | d_dC + y / dC\ + Vu> dC 
\ dt du \ du. du 



+ V-P* 



n 



v 'om + llenl1 + ^ + v ' r ' 



(18) 



where the tensor P* denotes the FLR-corrected CGL pressure tensor 



P* ee P CGL + eQ Pl _ (I - bb) = p|| bb + p ± (1 + efi) (I - bb), 



(19) 



the Hermitian conjugate operator 
field 

and the vector field 

T ee uA + ( P* 



II 1 " is defined in Eq. (3), and we introduced the scalar 
A . (20) 



dC 



+ 



Vj_ - — Vj 



r <9(V0) 
9(en||0||)' 



+ 



(21) 



9(Vi0) r 9(Vi0) 

Note that, while the last two terms d t A+ V ■ T in Eq. (18) do not play a role in the variational 
principle / 5C d 3 x dt = 0, they play a crucial role in the derivation of exact conservation laws 
based on the Noether method. 
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B. Gyrofluid Equations 



1. Gyrofluid velocity 



From the Eulerian variation of the Lagrangian density (18), the stationarity of the action 
functional with respect to a arbitrary virtual fluid displacement £ yields the Euler-Poincare 
equation 

d dC 







dt du 



+ V.^u-j + Vu.- + V-P* 



n 



+ eftVTj. 



(22) 



which describes the evolution of the gyrofluid velocity u. Upon substituting the Lagrangian 
derivatives (16) into the Euler-Poincare equation (22), and using the fact that the back- 
ground magnetic field B is assumed to be a time-independent nonuniform vector field, we 
obtain 

'dun 



= mn 



en 



+ u • Vum b uxB* + en Vib 

at c 



+ Vp± + T ± V(enfi) + V [(p A - efip ± )bb 
where p A = P\\ — P±, and we have introduced the following definitions 
B* = B + (rac/e) u\\ V X b 



(23) 



B* • b = B + (mc/e) My b- V X b >, 



b* = W/B\ = b + mu\\ C 



(24) 



with the magnetic-curvature term C defined as C = (c/eBjj) b X (b- Vb). 

Eq. (23) can also be written in a more compact form by introducing the gradient of ip 
evaluated at constant perpendicular temperature, V t t/>, defined through Eqs. (2) and (14) 



as 



l|V0|| + 



VB + 



VT, 



where d\\<j>\\/dB 



OB dT_ 
2 (T ± /B) tt. Then Eq. (23) can be written as: 



+ n vtl, 



= mn 



'du\\ 
~dt 



+ u- V«n b 



en 



uxB* + en V T i) + V-P*. 



(25) 



Note that the vector field B* is NOT a divergenceless field (since Vtty 7^ 0), as is common 
to all guiding-center and gyrocenter Hamiltonian models (in which parallel gyrofluid velocity 
u\\ is replaced with the parallel guiding-center velocity v\\). 

It is easy to see that Eq. (23) can be divided into two equations: one equation that 
expresses the gyrofluid velocity u in terms of the gyrofluid moments (n,u\\,p\\,p±) and the 
electrostatic potential 0, and one equation that describes the time evolution of the parallel 
gyrofluid velocity u\\. The first equation is obtained by taking the cross-product of Eq. (23) 
with b, which yields the following expression for the gyrofluid velocity 



U = Mil b* + 



cb 



en B 



- X < e n W> + Vp± + T± V (e n SI) + V (p A - efi p±) bb 



(26) 



Therefore, according to Eq. 26, the gyrofluid velocity u, consists of the parallel velocity 
■uii = u - b and the following perpendicular gyrofluid velocities 



u D = (c/enBS) b X Vp A 



= (c/5J)bxV^ 



(27) 



uc = (PA/n) C 

which represent the diamagnetic velocity, the generalized E x B velocity, and the curvature- 
drift velocity, respectively, and the following gyro-potential FLR corrections 

cT, b 



WS! = 



nBl 



X V(nfi) and w c = - ettT ± C 



(28) 



to the E x B velocity and curvature-drift velocity, respectively. 



2. Evolution equation for u\\ 



The evolution equation for parallel gyrofluid velocity u\\ can be derived by taking the 
dot-product of Eq. (23) with b*: 

( dun 



mn 



dt 



+ u-Vuii 



- b* 

- V 



en Vip + Vp ± + T ± V{enQ) 
(p A - ettp ± ) b 



(29) 
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where we have used the identity (valid for any function /) 

V-(/bb) = [v-(/b)]b + /(b-Vb), 

so that we find b* • V • (/ bb) = V • (/ b). The gyrofluid equation (29) for u\\ includes terms 
associated with the parallel electric field and its FLR corrections as well as parallel thermal 
forces. 

The set of gyrofluid equations of motion for the four gyrofluid moments (n,p\\,p±, u\\) are, 
thus, given by Eqs. (6), (7)- (8), and (29), respectively. Each of these gyrofluid equations 
involves the advection operator u • V and the divergence V • u; it is the diamagnetic part 
of the advection operator that must be eliminated from Eqs. (7)-(8) and (29) by adding 
suitable diamagnetic fluxes. 



C. Polarization Equation 

The polarization equation for <fi can be found by considering the terms in the variation 
(18) that involve the variation of the potential 5(p. Here, we must remember to add up the 
contributions of all the species to the total Lagrangian density. 

Using the variation of the Lagrangian density (18), we obtain the Euler-Lagrange equation 
for 0: 

= E IIM.f + V "9(V0)> (3°) 
where summation over particle species (J2j) is shown explicitly. Using Eq. (5), we find 

so that the Euler-Lagrange (30) becomes the polarization equation 

E 



+ -^V 2 = 0. (32) 



47T 



Comparing Eq. (32) with our previous one-temperature model [7], it is clear that the polar- 
ization equation is not affected by the introduction of the parallel temperature in the model. 
This arises from the fact that the terms that contain the gyrofluid electrostatic potential (2) 
do not depend on p\\. Further details concerning Eq. (32) can be found in Ref. [7]. 

Lastly, we note that two versions for the function V(x), appearing in the definition of 
the gyrofluid scalar potential (2), are commonly used in gyrofluid applications: the function 
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V(x) = e x or its Pade approximant V(x) = (1 — x)^ 1 . It is important to note that, 
whatever form is adopted for the gyrofluid scalar potential \\</>\\, it must be used consistently 
throughout the model in order to ensure energy conservation. 

IV. ENERGY CONSERVATION LAW 

In this Section, we present the local and global forms of the energy conservation law, 
as they arise from the application of the Noether method. For this purpose, we point out 
that, as a result of the variational principle / 5Cd 3 xdt = 0, the only remaining terms in the 
Eulerian variation of the Lagrangian density (18) become the Noether equation 



where A and T are defined in Eqs. (20)-(21). The energy and momentum conservation 
laws are derived from the Noether equation (33) by considering infinitesimal time and space 
translations, respectively. In the present work, we focus our attention on the local and global 
energy conservation laws associated with our electrostatic gyrofluid model. 

A. Local energy conservation law 

We derive the local form of the energy conservation law from the Noether equation (33) 
by considering infinitesimal time translations t — > t + St, from which we obtain the following 
expressions for the virtual fluid displacement £ and the Eulerian variations S(f> and SC: 



sc. f + v.r, 



(33) 



£ = — uSt, 5<p = —Std t 4>, and 5C = — St d t C. 



(34) 



Inserting these expressions into Eqs. (20)-(21), we obtain 





(35) 



where summation over fluid species is implied wherever appropriate. By combining these 
expressions, we arrive at the primitive form of the local energy conservation law; 



de' 

f + v.s^o, 
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where the primitive energy density is 

dC „ 1 



u • 



<9u 



2 P\\ ( , l E | 2 

L = -nmu\\ + pj_ + + [enip - — — 

Z Z \ 07T 



(37) 



and the primitive energy-density flux is 



+ 



V i ,x — V 



9(en||0||) 



+ 



dt ^ d{V\(f) x dt d{Vl4>) 
= u^nm-u 2 + p± + ^ + en-0^ + P* • u + 
d(en ||0||) \ 



dC 



+ 



dt V± V 9(Vi0) 



at <9(V0) 

90 / <9(en 



+ 



(38) 



9* V «9(Vi0) 

In order to arrive at the final form of the energy conservation law, we need to rearrange 
terms in Eq. (37). By substituting the polarization equation (30) into Eq. (4), we find 

IEI 2 



en 



where 



D = 



dC 



4:71 

d(en H0II) 



<9(V0) <9(Vi0) 
so that the last terms in Eq. (37) become 



+ ran Ug 



V± + 0V_ 



V-D, 

/ d(en 



9(V 2 



+ 



enij) 



IEI 



IEI 



8tt = 2 nm|Us| + Svr 



V-D. 



(39) 



(40) 



(41) 



Hence, we express the primitive energy density (37) as e' = e — V-D, where the final form 
of the energy density is defined as 

e = ^nm (u\ + |u E | 2 ) + p ± + ^ + 

and we obtain the local form of the energy conservation law 



(42) 



de „ ^ 



o, 



where the final form of the energy density flux is defined as 

After some partial cancelations, the final form of the energy density flux is 



uQnmw 2 + p± + y + entpj + P* • u + S^, 



(43) 



(44) 



(45) 
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where P* ■ u = p± (1 + eQ) u + (p& — eVLp^) u\\ b and the electrostatic energy density flux is 
defined as 

s --^fe-™ + ...) + v^|(^||i) + . ..).(«, 

In the next Section, after heat fluxes are inserted back into the pressure evolution equations 
and diamagnetic cancellations are performed with the addition of terms in the gyrofluid 
equations, the local energy conservation law (43) is converted into a new energy equation 
d t e + V ■ S* = 0, in which heat fluxes and diamagnetic-cancellation terms result in a modified 
energy density flux S*. This new form ensures that the total energy S = J ed 3 x satisfies the 
global energy conservation law dS/dt = 0. 



B. Global Energy Conservation Law 



The global energy conservation law d£/dt = can be derived from the local energy 
conservation law (43) by integrating it over space. Here, the global energy is defined as 



/ 



d 3 x 



p\\ 

nmu| + lp± + 



II IEI 



(47) 



In this form, the parallel kinetic energy, the internal energy, and the electric field energy 
explicitly appear. In a later section, we will present the time evolution of each of the separate 
terms that constitute the energy conservation law in order to identify the energy-exchange 
processes that allow the transfer of energy between the three types of gyrofluid (parallel 
kinetic, internal, and field) energies. 



V. DIAMAGNETIC CANCELATIONS AND ENERGY CONSERVATION 

The gyrofluid velocity (23) contains the diamagnetic velocity u D . Since the gyrofluid 
moment-equations are derived by inserting the gyrofluid velocity u in the Lagrangian con- 
straints, diamagnetic advection terms appear in the equations of evolution for the parallel 
velocity and the two anisotropic pressures. More specifically, the momentum equation for u\\ 
contains the term uo • Vizy, and the parallel and perpendicular pressure equations contain 
a combination of the terms up • Vp and pV • ud (where here p is either p\\ or p±). These 
diamagnetic-advection terms should be canceled in the gyrofluid evolution equations by the 
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introduction of appropriate terms containing higher-order moments [7]. These correspond 
to the FLR corrections to perpendicular fluxes in conventional fluid models [15, 16], but 
arise naturally due to grad-B and curvature drifts in the moment-based derivation of local 
gyrofluid models [4]. 

Since the higher-order moment terms cannot be derived from the Lagrangian, the dia- 
magnetic cancelations must be done manually. However, there exists a constraint in the 
addition of higher-order moment terms, namely, that the global energy conservation law 
should not be altered. Here, the terms added are derived from the Vlasov equation, and 
FLR corrections are then introduced to conserve energy. Thus, the final moment equations 
are not directly derived from the Lagrangian, but they still conserve energy exactly. 



A. Parallel Gyrofluid Dynamics 

The diamagnetic cancelation needed for the parallel momentum equation (29) involves 
the addition of the term — V-IIjj, associated with the non-diagonal part of the pressure 
tensor [13], on the right side of Eq. (29). From Vlasov theory, the diamagnetic-cancelation 
term is found to be [14] : 

P± X V«|| = mnu D • Vu|| + p± K,(mu\\), (48) 
eB \\ J 

where the magnetic differential operator /C(- • •) is defined by the identity 

V "(sJ* XV/ ) = ^-V/xVs - glC(f), (49) 

valid for arbitrary functions / and g. As written here, Eq. (48) conserves energy by itself, 
since u\\ V-IIjj = V- (u\\ Ily). Here, we should note that the substitution of B to B*\ is 
done here, and throughout this Section, to make the diamagnetic cancelation exact. 
The time evolution of the parallel kinetic energy density is, therefore, expressed as 

d firm 2 \ fmn , t-t*\ 

- u- [enV^p + Vp± + T ± V(enn) + v[(p A -eQp ± ) bb] } , (50) 

where we used the expression (26) for the gyrofluid velocity to write u\\ b* •{•••} = u ■{•• •} 
for the second term on the right side of Eq. (50). 
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B. Internal Energy 



To consider the time evolution of the internal energy p± + p\\/2, we rewrite the two 
pressure equations (7)-(8) in the form 

+ V-(p,|u) + 2 P ||bb: Vu = -2V.q,| ± + 4 q,^ • (b • Vb) + 2Q h (51) 
^ + V-(p ± u) + p ± (I-bb):Vu = - V-q ±± - 2 q,,^ • (b • Vb) + Q ± , (52) 

where Q\\ and are additional terms (to be determined later) associated with FLR cor- 
rections to the electrostatic scalar field 0. 

To zeroth order in the electrostatic potential 0, the heat fluxes are derived from the 
Vlasov equation directly and are found in Ref. [9] to be expressed as 

1 cp ± b 

qy± = 2 x VT II + H Uc ' ( 53 ) 



and 



q±± = 2 ^ X VT ± . (54) 



Inserting these diamagnetic heat fluxes, using the definition (49), we find 

- 2 V ■ q|| j. + 4q||j_-(b-Vb) = nu D • V7]| + p ± /C(7]|) 

- 2 p ± C • VTj, - 2 V • (t,| p A C) (55) 

for the parallel pressure equation (7), and 

-V-q ±± - 2 q NX - (b • Vb) = 2nu D -VT ± + 2 P± K(T ± ) +p ± C-VT ll (56) 

for the perpendicular pressure equation (8). 

Using Eqs. (51) and (52), the time evolution of the internal energy can, therefore, be 
expressed as 



^ + Rl) = - V • + p ± ) u + P CGL • u + (q,| ± + <u ± ) 

+ u-V-Pcgl + (Q|| + Q±), (57) 



where suitable energy-conserving expressions for the FLR-corrected heat fluxes Q\\ and Q± 
are now determined by considering the expression for the time evolution of the electrostatic 
field energy. 
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C. Electrostatic Field Energy 



By making use of the electrostatic field energy equation (41), we write the following 
expression for the time evolution of the electrostatic field energy 



d ( mn 



IEI 



dt 



2 + 8. 



„ «9D d ( , |Ep 



(58) 



where we made use of the polarization equation (32) and the definition (46) for S^. Here, 
using the gyrofluid continuity (6), we obtain 

on 

e (V> - T ± Q) — = - V • [ e {nil) - p± ft) u ] + enu ■ [ Vip - V(T ± ft) ] , 

and, using the perpendicular pressure equation (52), we obtain 
dp. 



dt 



= - eft V • (p ± u) - P n : Vu 



e n 



V-qi ± + 2q|| ± -(b-Vb) - Q_ 



= - V • ( en P± u + P n • u ) + u-[p± V(eft) + V ■ P n ] 



eft 



V-<u± + 2qn ± .(b-Vb) - Q_ 



where Pn = P* — Pcgl = eftp_i_ (I — bb) denotes the FLR-correction to the CGL pressure 
tensor. Hence, the first two terms on the right side of Eq. (58) can be written as 



- V • ( e nip u + Pn • u ) 
+ u-(enV^ - enftVT ± + V • Pn) 



e n 



V-q±± + 2 q||j_ • (b • Vb) - Q± . (59) 



We now require that, in order for the last term Qq = Q\\ + <3_i_ appearing on the right side 
of Eq. (57) to cancel the last three terms on the right side of Eq. (59), the latter terms must 
be written up to an exact spatial divergence as 



en V • q_L± - 2 en qin • (b • Vb) + eft Q_ 



V-Sc 



Hence, we set the additional terms Q\\ and Q± in Eqs. (57) and (59) to be 

Qi, = 2 eft qy ± - (b-Vb) 

Q± = -q±±-V(eft) - V- (q±±eft) 



(60) 
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so that 



Qn = Q\\ + Q± = nw c -V2|| + 2eQ P± IC{T ± ) + 2n {2w Q -eQu D ) -VT ± , (61) 

and = eVt (2 + efl) qj_j_, with q_|_j_ defined in Eq. (54). Combining Eqs. (59)-(61) into 
Eq. (58), the time evolution of the electrostatic field energy is, therefore, expressed as 

d (mn . l9 |E| 2 \ „ . , „ „ s 

— I — |u s | 2 + t-J-\ = - V- (en^u + P Q -u + + S n ) 

+ u- (enV^ - enfiVT ± + V • P n ) - Qn- (62) 



D. Explicit form of the energy conservation law 

When diamagnetic cancellations and heat fluxes are introduced into the four-moment 
gyrofluid equations, the local energy conservation law (43) is modified. By combining the 
evolution equations for the parallel kinetic energy (50), the internal energy (57), and the 
electrostatic field energy (62), the local energy conservation law (43) becomes the local 
energy equation 

de 



H + v-s = - V 

at 



u\\ m + q,| ± + (1 + en) 2 q ±± , (63) 



which ensures that the total energy £ = J ed 3 x still satisfies the global energy conservation 
law dS/dt = 0. 

We now identify the energy-exchange processes that transfer energy between the three 
different types of gyrofluid energy. First, we write down expressions describing the time 
evolution of each type of energy (e.g., parallel kinetic energy, internal energy, and field 
energy). Thus, the contribution of each species to the integrated parallel kinetic energy (50) 
is 

Jtl (iT^)^ = - / u ' ( enV ^ ~ enfiVT ± + V-P*)d 3 x, (64) 
the contribution of each species to the integrated internal energy (57) is 



d_ y/pji 
dtJ V2 



+ p ± ) d 3 x = J [ u • ( V • P CGL ) + Qn ] d 3 x, (65) 
and the contribution of each species to the integrated electrostatic field energy (62) is 

J ^ \u E \ 2 + ^Pjd 3 x = J[u> (enW - enfiVT ± + V • P n ) - Qn } d 3 x. 

(66) 

18 



d_ 

di 



The terms on the right side of Eqs. (64)- (66) appear in pairs with opposite sign, and give 
the energy-exchange processes. For example, the FLR-correction heat flux Qq is involved 
in energy exchange between the electrostatic field energy and the internal energy, while the 
FLR-correction pressure tensor is involved in energy exchange between the electrostatic 
field energy and the parallel kinetic energy. The contributions from en Vip, enQVT± and 
V • Pcgl in Eqs. (64)-(66), on the other hand, describe standard energy-exchange processes. 



E. Comparison with Previous Models 



We now write explicit final expressions for the gyrofluid density n, the gyrofluid parallel 
velocity u\\, and the parallel and perpendicular gyrofluid pressures p\\ and p±. The gyrofluid 
continuity is expressed in expanded form as 

dEii 



dt 



+ V 



n 



(u\\ b* + u c + w c )] = [/C(p±) + n/C(e^) + T ± lC(enQ)] 

+ u D -V(enfi) + wn-Vn, (67) 



where ds/dt = d/dt + u^, • V. With the insertion of the diamagnetic-cancellation term (48), 
the evolution equation for the gyrofluid parallel velocity is expressed in expanded form as 

/ dm 



mn 1-7^ + u'-Vuiij - p±K.(mu\\) = - b* • [ en + Vp± + T± V(enfi) 



- V 



(p A - eQp±) b 



(68) 



where u' = u — denotes the gyrofluid velocity without its diamagnetic contribution. 
Lastly, with the insertion of the diamagnetic-cancellation terms (55)-(56) and (60), the 
gyrofluid parallel and perpendicular pressure equations are 



^ = - V- (p|| u + 2p||«u b + 2q,| ± ) + 2u-V- (p,|bb) 
+ 4(l + efi)q||j..(b.Vb), 



(69) 



and 



dt 



- V' 



P± u + pj_ (I - bb) • u + (1 + eQ) q ±± 



P± (I-bb) 



- 2q|| ± .(b-Vb) - q±±-V(efi), 



(70) 



where the heat fluxes qin and qj_^ are defined in Eqs. (53) and (54). 
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The results of the present two-temperature gyrofluid model can be compared with the 
previous one-temperature model presented in Ref. [7], where the gyrofluid equations are 
derived by including the perpendicular temperature only. Since the FLR-corrected CGL 
pressure tensor (19) includes the parallel pressure p», the two gyrofluid models agree except 
for the terms that arise from the pressure anisotropy p& = p» — p±. 

Our results can also be compared with the gyrofluid model of Beer and Hammett [4] , since 
this Beer-Hammett model is the most extended one, which includes closures and contains 
all the previously developed gyrofluid models. The equations of evolution of the gyrofluid 
moments presented here are nearly identical to those of Beer and Hammett. The differences 
that arise can be separated in two categories: (1) differences in the non-FLR terms and (2) 
differences in the FLR terms. 

The non-FLR terms coincide almost exactly in the two models, except terms that come 
from closures of higher-order moments, which our model is not able to retrieve. In our 
model, gyrofluid moments higher than pressure moments cannot be included in the gyrofluid 
Lagrangian, and the closures are done automatically when choosing the constraints for the 
variation. Thus, for instance, Landau damping is not included in our model, though it can 
be added by hand afterwards as long as the energy conservation remains exact. A special 
difference in the non-FLR terms between the two models - and the only one of the kind 
that appears - is a different term in the momentum equation (68): the magnetic term JC(u\\) 
in our model compared to 2JC(u\\) in Ref. [4]. This difference arises from the fact that the 
additional magnetic contribution originates from the parallel-parallel heat flux qyn, which 
we cannot fully retrieve. 

The differences in the FLR terms between the two models arise from the higher-order 
moment closures, but also because of the energy conservation law. Especially in the FLR 
terms, previous models do not conserve energy. Although our polarization equation reduces 
to the local one [5] and so is in agreement with the local model in the local limit, we do 
not exactly retrieve the results of Ref. [5]. That is, although we find a strong correlation 
between n and ||0|| and p± and Q, there exists another FLR term that is not correlated with 
p±. This difference arises because of the constraints introduced here, which lead to another 
closure than that of Ref. [5] . 
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VI. SUMMARY AND CONCLUSIONS 



In this paper we have derived a set of electrostatic gyrofluid equations for an anisotropic 
plasma, that describe the evolution of density, momentum, parallel and perpendicular pres- 
sures. The fully inhomogeneous four-moment model also includes a polarization equation, 
from which the electric field is computed, and satisfies an exact energy conservation law that 
includes energy-exchange terms involving parallel kinetic energy, internal energy, and field 
energy. 

To guarantee energy conservation, we have used the Lagrangian approach in which all the 
equations are derived by a variational principle. Diamagnetic cancelations were taken into 
account in a second step, by the addition of terms that include gyrofluid moments higher 
than the pressure. All non-electrostatic diamagnetic-corrections terms were computed from 
the Vlasov equation, and FLR-corrected terms associated with <p were chosen so that the 
energy conservation law is still satisfied. 

The validity of the Lagrangian and the constraints introduced for the variational proce- 
dure was verified by the validity of the final evolution equations and the energy conservation 
theorem. The fact that the set of equations conserves energy makes the model suited for 
turbulence computations, and especially for those that treat large-amplitude disturbances, 
or those that treat strong variations of the plasma parameters. Moreover, since the set 
of equations is fully inhomogeneous, it is suitable for studying the nonlinear evolution of 
various fields, independent of the linear growth phase. 

The successful application of this technique to the four-moment gyrofluid model show 
that it is feasible, although not trivial, to extend the Lagrangian formulation to include also 
the heat fluxes as dynamical variables, and, thus, arrive at a complete six-moment model. 
This extension will be treated in future publications. 
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